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Abstract. In the case of a one-dimensional nonsingular Hamiltonian H and a singular 
supersymmetric partner H a , the Darboux and factorization relations of supersymmetric 
quantum mechanics can be only formal relations. It was shown how we can construct 
an adequate partner by using infinite barriers placed where are located the singularities 
on the real axis and recover isospectrality. This method was applied to superpartners of 
the harmonic oscillator with one singularity. In this paper, we apply this method to the 
singular isotonic oscillator with two singularities on the real axis. We also applied these 
results to four 2D superintegrable systems with second and third-order integrals of motion 
obtained by Gravel for which polynomial algebras approach does not allow to obtain the 
energy spectrum of square integrable wavefunctions. We obtain solutions involving parabolic 
cylinder functions. 
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1 Introduction 

In recent years, many papers [2, 4, 5, 8, 13, 17, 19, 23, 38] were devoted to a nonsingular isotonic 
oscillator and its various generalizations. It is also referred in literature as CPRS system and 
also often written using the following parameter ui = ^ or taking h = 1. This quantum system 
is exactly solvable and given by the following equation (i.e. with the parameter given by a = ian, 
a G M) 



Let us also define the singular isotonic oscillator that correspond to the Hamiltonian given by 
the equation (1) with alternatively the parameter defined as a G K and the singularities are 
on the real axis. Let us mention that in the nonsingular case the solutions in terms of series 
was obtained [8] that involve polynomials that are linear combinations of Hermite polynomials. 
Moreover, it was shown by J. Sesma [38] that the corresponding Schrodinger equation can be 
transformed into the confluent Heun equation [39, 40] and studied quasi-exactly solvable analogs. 
In addition, the CPRS system was related to position-dependent effective mass Schrodinger 
equations (PDMS) [23] and it was shown recently, using Bethe ansatz method [2], that this 
system possesses also a hidden sl(2) algebraic structure. 

Let us mention that the nonsingular isotonic oscillator was obtained and studied in earlier 
works of Spiridonov and Veselov in the context of the dressing chain method [41, 43] and that 
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some results presented in [13] were obtained in the papers [37, 42]. This system is also a parti- 
cular case of a system [3, 6, 7, 26, 28] involving the fourth Painleve transcendent [20] and related 
to third-order supersymmetric quantum mechanics (SUSYQM). Recently, using a modified fac- 
torization a family that include this system was also obtained [13]. 

The nonsingular isotonic oscillator also appears in context of 2D quantum superintegrable 
systems with a second and a third-order integrals of motion with separation of variables in Carte- 
sian coordinates classified by S. Gravel [18]. Four of these systems that possess more integrals 
than degrees of freedom, contain the nonsingular isotonic oscillator. These four superintegrable 
systems were studied [27, 30] from the point of view of SUSYQM [16, 22, 45] and polynomial 
algebras. The square integrable wavefunctions and the energy spectrum were obtained. Fur- 
thermore, the superintegrability property does not depend on the nature of the singularity and 
thus the 2D superintegrable systems remain superintegrable when we take a£R and thus these 
2D superintegrable systems also admit two types, i.e. one singular and one regular type. The 
singular version of these 2D superintegrable Hamiltonian remain to be solved. 

Let us mention, that the singular isotonic oscillator (i.e. equation (1) with a G R) belongs 
to the class of ID Hamiltonian constructed by Robnik [36] from SUSYQM with higher excited 
states to construct the superpotential. It was recognized that the factorisation and interwining 
(or Darboux) relations of SUSYQM are only formal. 

It was observed by many authors that for systems with singularities on the real axis related 
to a regular superpartner, SUSYQM does not allow [9, 10, 21, 31, 46, 47] to relate the wave- 
functions and the energy spectrum of the superpartners using supercharges. Supersymmetry 
and it algebraic equations (factorization, Darboux or intertwining) are formal and do not take 
into account boundaries or singularities. A method to reobtain the isospectrality property, i.e. 
relate the wavefunctions and the energy spectrum for a singular Hamiltonian and its regular 
superpartner, consists to use infinite barriers to modify the regular Hamiltonian [31]. However, 
this idea was only applied to systems with only one singularity on the real axis. The study of 
Hamiltonians with many singularities using this approach is an unexplored subject. 

The purpose of this paper is to apply this approach to the singular isotonic oscillator which 
has 2 singularities of the real axis and also use these results to obtain the energy spectrum 
and wavefunctions for four 2D singular superintegrable systems with a second and third-order 
integrals that remains to be solved. 

Let us present the organization of the paper. In Section 2, we recall results obtained for the 
constrained harmonic oscillator. In Section 3, present definition and results concerning super- 
symmetric quantum mechanics and discussed how we can construct the adequate superpartner 
for singular and nonsingular superpartners. In Section 4, we apply the results of Section 2 and 
the method of Marquez, Negro and Nieto [31] to the singular isotonic oscillator. We obtain the 
wavefunctions and square integrable wavefunctions. In Section 5, we discuss how the polynomial 
algebra obtained for four 2D singular superintegrable systems with a second and a third integrals 
of motion that contain the singular isotonic oscillator is only formal and related to SUSYQM. 
We apply the results of Section 4 to solve these systems and obtain the corresponding energy 
spectrum and wavefunctions in terms of the parabolic cylinder functions. 

2 Constrained harmonic oscillator 

Let us recall known results concerning the ID unconstrained and constrained harmonic oscillator 
(i.e. with one or two symmetric infinite barriers). The unconstrained harmonic oscillator that 
we denote by V\ is given by the following well known Hamiltonian 
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— OO < X < oo. 



(2) 
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The energy spectrum and the wavefunctions are also well known 



where h n are the Hermite polynomials. 

We can defined on [B, oo) and [-B, B] respectively the constrained harmonic oscillators Vn 
and V m [11, 32]: 

2 2 

li^-^X' X ~ B ' 5> °' (3) 
x < B, Bel, 

Hn(x)=<-r- ~ B ^± B - B >°- (4) 
x<-B, x>B, Bel. 

We consider also the following transformation 

and the Schrodinger equation corresponding to the non vanishing part of these potentials Vn 
and Vin thus become the Weber differential equation 





d?y(z) ( z 



•2 



+ e)y(z) = 0. (5) 



dz 2 \ 4 

2.1 Constrained harmonic oscillator Vn 

For the quantum system Vn given by equation (3) we have the following constraints on the 
solution y(z) of the equation (5) [11, 31, 32]: 

(i) y(z) must be square integrable, 

(ii) y(z) must be continuous on K, i.e. y{z) — > when x — > b + . 

The equation (5) can be solved in terms of the parabolic cylinder functions [1, 44] 

/ \ _1 7 2 /l 1 1 1 ,\ . , 1 2 /l 3 3 1 2 

W (c,z) = e-i« i^i(2 e+ 4'2'2 Z J' V2 ^ z) = ^ ^U' + i'^* 

where \F\ is the confluent hypergeometric function also related to the Whittaker function. We 
can also introduce the following and more appropriate functions in the case of this constrained 
harmonic oscillator [1, 44] U(e, z) and V(e,z) that are combinations of yi(e,z) and y2(e,z) and 
given by 

U(e, z) = Y 1 cos (vr Q + ) " Y * sin (* (\ + \ e 
V(e, z) = — x (yi sin (n (\ + \e J J + Y" 2 cos f vr (i + \e 

where 



T(i-e) V V \4 2 J J V \4 2 



7T2 sec (it (\ + ie)) , . 7T2 cosec (vr + |e)) 

2rtr(f + ie) 2 5 £ -irQ + |e) 
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The asymptotic behavior of the functions U (e, z) and V(e, z) is given by the following expressions 



,2 



V(e, x) diverge for z — > oo but U (e, x) is an appropriate solution and we thus take 

ju(e,z), z>b, 
[0, z < b. 

Let us present some results [11] concerning the discrete ensemble of solutions. It can be 
shown (from U(e%(b),b) = 0) that the energy levels depend of the position denoted b of the 
infinite barrier ( with b assumed to be positive ) in the following way 

enO) = -g + £ o( n ) + e i( n ) 6 + £ 2(n)fr 2 + • • • , n = 0, 1, 2, . . . , 

and the first coefficients are given by the following recurrence relations 

eo(n) = -(2n + 1), ra = 0,1,2,..., 

2n + 3 . . . , i i 

£l(n + 1) = 2^T2 £l(n) ' ei(0) = -25vr-2, 

r j.n /2n + 3\ 2 , , , (2n + 3)!(2n + 2)! 
e 2 (n + 1) = - — —r e 2 (n) 



2n + 2j w 16vr((n + l)!) 4 (n + 1)2 4 "' 
e2 (0) = -2^- 1 (l-Log[2]). 

These recurrence relations can be solved to obtain explicitly the value of the coefficients. The 
coefficients of higher-order terms can be calculated and this calculation could also be imple- 
mented numerically. Therefore, for every fixed value of b where the infinite barrier is located, the 
zeros of this equation provide a discrete ensemble of solutions [11, 31, 32] eJJ(i>), n = 0, 1, 2, . . . 
and for the Hamiltonian with the potential given as in equation (3) we have E^(b) = —hu)€n(b). 
The corresponding eigenfunctions are then U(e^(b),z). The Fig. 1 presents the energy for the 
ground state and the 10 first excited states as a function of the position of the infinite barrier. 
The Fig. 3 presents the ground state for an infinite barrier located at b = 1. The only exactly 
solvable case [11, 31, 32] is when 6 = 0. In this particular case we obtain from the condition 
U(e I l0) = 



E„ =-tkje n = fkj ^2n+l + , n = 0,1,2, 



The wavefunctions correspond to the odd levels of the harmonic oscillator. The limiting case 
where b — > oo allow to recover the spectrum of the harmonic oscillator given by equation (2). 

2.2 Constrained harmonic oscillator Vm 

For the system V\u given by equation (4), we have the following constraints [11, 31] 
(i) y(z) must be square integrable, 

(ii) y(z) must be continuous on M, i.e. y(z) — > when z — > b + and z — > —b~. 
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The associated eigenfunctions are the even and odd functions y±(e, z) and Jfefe z). 
2/m(e,z) : 



yi(e,z), n odd, —b<z<b, 
yiif-, z), n even, —b<z<b, 
0, z < -b, z > b. 



The expansion of the eigenvalues as a function of the positions of the barriers takes the form 
in this case [11] 

en'O) = + eo(n) + e 2 (n)b 2 + e 4 (n)6 4 + e 6 (n)6 6 + • • • , n = 0, 1, 2, . . . , 

with the first coefficients given by 
n + 1\ 2 



e„ 2 (n) = -(^— j tt z , n = 0,1,2,... 

e (n) = e 4 = 0, e 2 (n) = — — - . ■ 

12 8e_2(n) 

1 5 7 

e e( n ) = ^t; — t r + tf^. — r~v? + 



720e_ 2 (n) 192e_ 2 (n) 2 128e_ 2 (ra) 3 ' 

These relations thus provide for a fixed value of b a discrete ensemble of solutions eJJ I (6), n = 
0,1,2,.... For the system given by the equation (4) the energy spectrum is thus E^ l (b) = 
—hujc~ (b). The corresponding eigenfunctions are then yi(e„ (£>), z) for the even levels and 
yx(e~ r(6), z) for the odd levels. The Fig. 2 presents the energy for the ground state and the 10 
first excited states as a function of the position of the infinite barrier. The Fig. 4 presents the 
ground state for an infinite barrier located at b = 1. 

3 Factorization method and singular Hamiltonian 

In supersymmetric quantum mechanics [16, 22, 45], we can define two first-order operators A 
and (supercharges) in the following way 

where W(x) is called the superpotential. 
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Figure 3. The wavefunction for the ground 
state U(e^,z) (with h — 1 and b — 1) for the 
Case Vij. 




Figure 4. The wavefunction for the ground 
state 2/1(6™, z) (with h = 1 and b = 1) for the 
Case Vin- 



We consider the following two Hamiltonians which are called "superpartners" 

It can be shown that these algebraic relations allow to relate eigenfunctions and energy spectrum 

of the two superpartners H\ and H 2 . There are two cases that can be considered. The first one 

corresponds to broken supersymmetry and happen when Aij)\p 7^ 0, E^ 7^ 0, A^ipQ 7^ and 
(2) 

Bq / 0. The wavefunctions and the energy spectra are thus related in the following way 

Eg* = EW > 0, </4 2) = (E^A^, ^ = (JS^)-MVW, ( 8 ) 

and the two Hamiltonians are strictly isospectral. 

For the second case the supersymmetry is unbroken and one of the Hamiltonians allows 
a zero mode, i.e. a states is annihilated by one of the supercharges. Without lost of generality 
we take H\ as having a zero mode (also referred as zero energy ground state) and thus we have 
Aip^ = 0, Eq 1 ^ = 0, A^ip^ 7^ and E^ 7^ 0. We obtain the following relations between the 
wavefunctions and the energy spectra 

= Ei% EM = 0, ^) = (E^A^ n % tfi = (4 2) r^ 2) . 0) 

In this case, the superpotential W(x) can be written in term as W(x) = ~(iy- 

These relations (8) and (9) allow, if one of these systems (H\ or H2) is exactly solvable, to 
obtain the energy spectrum and the wavefunctions of its superpartner. A such construction can 
be extended more generally to families of superpartners and superpartners of £;th-order. It can 
be shown, that in the case of Hermitian operators, the factorization relations (7) are equivalent 
to the following intertwinning (also called Darboux) relations 

H X A = AH 2 , A ] H X = H 2 A\ (10) 

The factorization and interwining relations can be combined and written in terms of matrices 
that generate a superalgebra. These algebraic relations also allow to relate the creation and 
annihilation operators of the Hamiltonians H\ and H 2 ■ The ladder operators and c of H\ and 
the ladder operators of H 2 that we denote and M are related by 

M = A^cA, M ] = A^c^A, (11) 
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However, as mentioned earlier the algebraic relations given by equations (7), (10) and (11) 
can be formal and do not take into account boundaries and singularities. In the case of a ID 
nonsingular Hamiltonian H\ and a singular supersymmetric partner H2 the wavefunctions ip( 2 > 
obtained from those of H\ (i.e. ij)^) and by acting with the supercharges (i.e. using (8) or (9)) 
on are not guaranteed to belong in the Hilbert space of Hi- If the supercharges admit singularity 
the resulting wavefunctions could also admit singularity. 

It is known that a singularity of the form r^z^p makes that a particle is confined in one of 
the regions (i.e. left or right side of the singularity) because the probability of going through 
the barrier is zero. We thus need to compute independently the eigenfunctions for each region 
divided by a such singularity [15, 24]. It was shown by Marquez, Nieto and Negro [31] that 
for Hamiltonians with such singular terms we can construct an adequate superpartner by using 
infinite barriers placed where are located the singularities on the real axis and recover isospec- 
trality between the superpartners. The wavefunctions and the energy spectrum of the systems 
with singularities can now be obtained from those of the regular Hamiltonian with infinite bar- 
riers and the equations (8) and (9). Using this approach, families of quantum systems that are 
superpartners of the harmonic oscillator with one singularity were studied [31]. 



4 Singular nonlinear oscillator 

Let us consider the following supercharges that admit singularities at a and —a 

\/2 \ dx 2a \x — a x + a 
h f d 1 ( -1 -1 



A = —= —rx - + 



y/2 \dx 2a 2 \x — a x + a / 
We can form from these operators, the following Hamiltonian H\ and H2 

*, = ^ = f +w = £ + |£_l* m 

These algebraic relations can be used in the case a = iao, ao £ K to solve the nonsingular 
isotonic oscillator, but in the singular case o£i they can not be used. The Hamiltonian H2 
is the well known harmonic oscillator with an additive constant and in principle the SUSYQM 
relations would provide the corresponding wavefunctions and energy spectrum of the Hamil- 
tonian H\. However, the operators contains singularities and using SUSYQM relations, we 
obtain eigenstates of the Hamiltonian H\ but non square integrable states and thus non physical 
states. However, as mentioned this is possible to recover the isospectrality by using an appro- 
priate superpartner for each region we would like to compute the wavefunctions and defined by 
the singularities, i.e. R\ = (—00,— a], R2 = [a, 00) and R3 = [—a, a]. Due to symmetry, the 
calculation for the first and the second regions are similar. Let us only present, the calculation 
in the case of the regions R2 and R3. 

In the case of the singular Hamiltonian H\ in the region R2 we consider a modified super- 
partner that is the potential given by equation (3) up to an additive constant that does not 
affect the results obtained in Section 2 is 

h 2 x 2 5h 2 

V 2u (x)= { ~8^~4^' a> °' 

x < a, a£l. 
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In the case of the singular Hamiltonian Hi in the region R% we take a modified superpartner 
as given by equation (3) up to an additive constant 



V2 1U (X) 




—a < x < a, a > 0, 
x < —a, x > a, a £ 



The case with two singularities differs from the case with one singularity because of the 
possibility of bounded intervals as i?3 and as we saw from the Fig. 1-4 the behavior of the 
two types of constrained harmonic oscillator differ and thus their superpartners will also admit 
different spectrum and wavefunctions. We can now apply results of Sections 2 and 4 and we 
obtain for the Hamiltonian given by equation (12). 

Case i?2 = [a, oo): 



rP(e^{a),x) = A^y n (x) 
H / 



d 



_y/2 V dx 2a 2 ' 



-1 -1 

+ 



x — a x + a 



U{e l n{a),x), x>a, 
x < a, 



E, 



h 2 

2„^ f;;U 



Case -R3 = [—a, a]: 



E 



in 



a),x) 

h_ 
h_ 
o, 

~ 2a 2 



= A ] ym(x) 
d 1 
dx 2a 2 
d 1 
dx 2a 2 



-1 -1 

+ 



x — a x + a 

-1 -1 

+ 



x — a x + a 



yi (e„ (a), x) , n even, x < \a\, 

y 2 (e 1 r l 1 (a),x), n odd, x < \a\, 
x > lal, 



4 n («) 



The wavefunctions are square integrable and also continuous at the infinite barrier by construc- 
tion. This can be obtained by considering the local behavior in terms of a Taylor series about 
the position of the infinite barrier. The completeness of the set is also provided by the fact that 
the state annihilated by the operator A has singularities at the infinite barrier and is not an 
admissible wavefunction. In the next section we will show how these results can also be used in 
context of superintegrable systems that also involve such singular systems. 



5 SUSYQM and superintegrability 

Let us now consider the following 2D superintegrable Hamiltonian 

= H « + H « = + + fi2 (n** + j^w + raO ■ (14) 

This system is one of the four irreducible superintegrable systems with a second and a third- 
order integrals of motion found by Gravel [18] that remained to be solved. An important 
property of superintegrable systems with higher-order integrals of motion is that the quantum 
systems and their classical analog do not necessarily coincide [18]. In this case the classical 
limit is the free particle. It involves the nonsingular or the singular isotonic oscillator as we 
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choose a = iao, do £ I or a 6 1. The existence of the integrals of motion does not depend on 
the nature of the singularities and thus of this choice. This system was studied [27, 30] using 
polynomial algebras for both cases. However, only in the nonsingular case the finite-dimensional 
unitary representations correspond to physical states. In the singular case, the integrals (i.e. 
they commute with the Hamiltonian H s \) I\ and I 2 take the form 



" ' J ' " 1 8a 4 (si -a) 2 ( Xl + a) 2 

1 2 f f (x 2 -4a 2 ) 2 4(x 2 + a 2 ) \ 

+ 2 n T 1 V 4a 4 x \ - a 2 + (x? - a2)2 J 

where { , } is the anticommutator. They are respectively polynomials in the momenta of order 2 
and 3. These integrals generate the following commutation relations [27, 30] 



4/i 4 

[h,I 2 ]=I 3 , [h,h] = —I 2 , 

ar 



h 4 



[I 2 , h] = -2h 2 l\ - m z i'(H sl + 8^H^ + 6^1? 

+ 8—^H s iIi — 8—^H sl + 2— jJi — 2^-Hsi — 6-77. 
a z a z a 4 a 4 a° 

We can calculate the Casimir operator, the realizations in terms of deformed oscillator al- 
gebras, the finite-dimensional unitary representation and the corresponding energy spectrum. 
For a € R the following solutions is obtained 

E = HN) = - 4 N(p + 1 - N)(N + 1)(N + 3), 

where $(iV) is the structure functions of the deformed oscillator algebra. However, the polyno- 
mial algebra does not take into account boundaries and singularities and is thus only a formal 
algebraic construction in this case and this finite-dimensional unitary representations do not 
correspond to physical states. Furthermore, this is the same energy spectrum that is obtained 
by using the formal factorization and Darboux relations with the harmonic oscillator and given 
by equation (12) and (13). However, as mentioned we can show they are not square integrable. 

The fact that the polynomial algebra and the SUSYQM relation are formal is also related. 
This can be understood from the fact that the integrals and the polynomial algebras can be 
derived using SUSYQM [26]. In the xi-axis, the ladder operators are constructed from super- 
symmetry with the supercharges using equation (11) [26, 27]. The annihilation operators are 
given by 

h 2 ( d 1 (I 1 
M x 1 = — [ —, + ( + 



4a 2 \ dx\ 2a 2 \x\ — a x\ + a 



x I xi + 2a 2 — II- —=xi + 



d \ f d 1 / 1 1 



dx\ I \dx\ 2a 2 \x\ — a x\ + a 



h ( 2 d 



M - = 2^ r + 2 ° dx 2 

with the creation operators obtained by considering {M xi y and {M X2 y . The integrals can be 
written as 

— 2a 2 i —2a 2 i 
h = -j—iMlM^-M^Ml), I 3 = —^—{MlM^+M^Ml). 
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5.1 Modified SUSYQM and singular superintegrable systems 



Let us consider again the Hamiltonian H s \ given by the equation (14) and use the results 
of Section 4 on the ID Hamiltonian H Xl and takes a modified superpartner that is isospectral 
using infinite barrier at the singularities. The ID Hamiltonian H X2 is only the ID unconstrained 
harmonic oscillator. 

The wavefunctions and the energy spectrum would be thus for a particle in the region [a, oo) 
for the axis X\ 



E 



h 2 



h ( d 1 
y/2 V dx\ 2a 2 



Xl 



-1 



■1 



x\ — a x\ + a 



2a2 (£(«) + *" 2). 
For a particle in the region [—a, a] for the axis x\ 



$(xi,x 2 ) = x(x2)A^y^j(xi) 



X{X2) 
X{X2) 



h / d i 

y/2 V dxi 2a 2 



n 



d 



y/2 V dxi 2a 2 



xi 



xi 



-1 



■1 



xi — a xi + a 



-1 



xi — a 



+ 



-l 



xi+a 



E 



0, 

h 2 



2a 2 v n 



All, 



+ k-2), 



with 



x (x 2 ) = Ce^H k 



2a 2 



X2 



U(e n (a),xi), xi > a, 
x\ < a, 



yi(e n (a),xi), n even, x\ < \a\ 

y 2 (e n (a),^i), n odd, x\ < \a\, 
xi > lal, 



The eigenfunctions and energy spectrum for the three other irreducible 2D superintegrable 
systems with a second and third-order integrals of motion can be calculated using the same 
approach. 

Hamiltonian H S 2 



1 



1 



Hs2 = -P^ + -P Z £2 +^ 



2 / x\ + 9x?, 



+ 



1 



+ 



8a 4 (xi - a) 2 (xi + a) 2 
For a particle in the region [a, oo) for the axis xi 



<5>(xi,x 2 ) = X (x2)A^y£ ) (xi) 



E 



X{X2) 



h 2 



d 



l 



y/2 V ^i 2a 2 



-1 



1 



+ 

xi — a xi + a 



U(e n (a),xi), xi > a, 
xi < a, 



2a2 (e^(a) + fc-l). 
For a particle in the region [—a, a] for the axis xi 
$(xi,x 2 ) = x(x 2 )^ t 2/ I ( I a I ) (xi) 
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X{X2) 



-( 


' d 






v dx\ 




-( 


f d 




V2\ 


v dx\ 



-1 -1 
+ 



x\ — a x\ + a 



-1 



+ 



■1 



x\ — a x\ + a 



0, 

^=^( £ n I ( a ) + fc - 1 )' 



yi(e„(a),xi), n even, x\ < \a\, 

y 2 (e„(a),xi), n odd, xi < |a|, 
x\ > \a\, 



with 



x (x 2 ) = Ce^ x ^H k 



Hamiltonian 



2a 2 



X 2 



8a 4 x\ ^ {x\ — a) 2 ' {x\ + a) 2 



+ 



For a particle in the region [a, oo) for the axis x\ 



Hxi,x 2 ) = X (x2)^y^\x 1 ) 



X{x 2 ) 



h 



d 



V dxi 2a? 



xi - 



-1 -1 

+ 



x\ — a x\ + a 



0, 



?7(e„(a),xi), xi > a, 
xi < a, 



^=2rf (£(*) + *)■ 
For a particle in the region [—a, a] for the axis xi 



$(xi,x 2 ) = x(x 2 )^ t y I ( n ) (xi) 



X{X2) 
X{X2) 

0, 

h 2 



h 



d 



l 



v/2 V 2a2 



ft 



V dxi 2a 2 



xi 



xi - 



-1 



1 



xi — a ^ x\ + a 

-1 -1 

+ 



xi — a xi + a 



yi(e„(a), xi), n even, xi < |a| 

y 2 (e„(a),xi), n odd, xi < |a|, 
xi > |a|, 



* ^ (e " I(a) + /C )' 



with 



x (x 2 ) = e ^ x,i x\L\ 



2a 2 ' 



where are Laguerre polynomials. 
Hamiltonian H s ^ 



Hs4 = \PI + \P 2 X2 + fr 2 



X? + Xn 1 1 

+ 7 w> + 



r) 



1 1 

+ 



8a 4 ' (xi-a) 2 (xi + a) 2; (x 2 - a) 2 (x 2 + a) 2 J ' 
For a particle in the region [a, oo) for the axis xi and x 2 



*(x u x 2 ) = f[4yg>(x i ), 
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with 




d 1 



y/2 V dxi 2a 2 ' 



1 -1 

+ 



Xi — a Xi + a 



For a particle in the region [—a, a] for the axes x\ and x 2 

2 

t.,(«)/ 



C/(e n (a),Xj), > a, 
Xj < a, 



*(x 1) x 2 ) = nA t Kra(^). 



with 



' r h_ f d i_ 

7! \dx~ ~ 2a? Xi 

" h_ f d i_ 

7! \dx~ ~ 2a^ Xi 



+ 



Xj — a 



-1 -1 
+ 

Xi — a Xi + a 



0, 

h 

2a 2 



2/i(e n (a),Xj), n even, Xj < |a| 

y 2 (e n (a),2;j), n odd, < |a|, 
xi > lol, 



-/^^(d H (a) + ^-5), 



6 Conclusion 

In this paper, we studied the singular isotonic oscillator using the approach introduced in [31]. 
We obtained the energy spectrum and the wavefunctions from supersymmetry by identifying 
the appropriate superpartner. We also studied four singular 2D superintegrable Hamiltonians 
obtained by Gravel. The energy spectrum in the ID cases is no longer equidistant and in the 2D 
superintegrable versions the energy spectra do not display accidental degeneracy explained by 
a polynomial algebra that is only formal. However, this point out how supersymmetric quantum 
mechanics can be used to solve such systems. The wavefunctions are also interesting from the 
point of view of special functions as they are written in terms of expressions involving parabolic 
cylinder functions. These results could also be interesting in regard of possible applications. The 
constrained harmonic oscillators has found applications in the context of condensed matter [14]. 

Moreover, the factorization is not unique [33] and it was shown that the nonsingular isotonic 
oscillator has more general families of superpartners [29]. Many families of nonsingular super- 
partners of the radial oscillator were obtained and the relations with exceptional orthogonal 
polynomial studied [35]. The study of the singular version of these systems and the relations 
with special functions is also important. The approach described in this paper could also be ap- 
plied to systems involving the fourth and the fifth Painleve transcendents that generate families 
of systems with singularities on the real axis for specific parameters [3, 6, 7, 25, 26, 28]. 

These results are also important as supersymmetry plays a role in the construction of new 
superintegrable systems in Cartesian [25, 29] but also in polar coordinates [12, 34]. 
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